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Abstract 

We study the problem of covering a given set of n points in a high, d-dimensional space by 
the minimum enclosing polytope of a given arbitrary shape. We present algorithms that work for 
a large family of shapes, provided either only translations and no rotations are allowed, or only 
rotation about a hxed point is allowed; that is, one is allowed to only scale and translate a given 
shape, or scale and rotate the shape around a fixed point. Our algorithms start with a polytope 
guessed to be of optimal size and iteratively moves it based on a greedy principle: simply move 
the current polytope directly towards any outside point till it touches the surface. For computing 
the minimum enclosing ball, this gives a simple greedy algorithm with running time 0(nd/e) 
producing a ball of radius 1 + e times the optimal. This simple principle generalizes to arbitrary 
convex shape when only translations are allowed, requiring at most 0(l/e 2 ) iterations. Our 
algorithm implies that core-sets of size 0(1/ e 2 ) exist not only for minimum enclosing ball but 
also for any convex shape with a fixed orientation. A Core-Set is a small subset of poly(l/e) 
points whose minimum enclosing polytope is almost as large as that of the original points. 
When only rotation about a fixed point is allowed, for a certain class of convex bodies with 
an axis of symmetry that includes cylinders, cones and ellipsoids, we prove that our techniques 
work provided the problem is confined to a half space. Without the half-space restriction, we 
obtain an algorithm whose running time is exponential in 1/e 2 , and corresponding core-sets. 
This automatically gives us an 2°( 1 / e ^nd time algorithm for the min-cylinder problem provided 
we are given a fixed point on the axis. Although we are unable to combine our techniques for 
translations and rotations for general shapes, for the min-cylinder problem, we give an algorithm 
similar to the one in but with an improved running time of 2°^ los *'nd. This generalizes 
to computing the minimum radius fc-dimensional flat in time exp( e -g - log ^)nd. 

1 Introduction 



Given a set S of n points in d dimensions, we study the problem of finding the minimum enclosing 
polytope of a given arbitrary shape when d is large. Being a fundamental problem in computational 
geometry with applications in data mining, learning, statistics and clustering ([5], [5], |S]), this 
problem has a rich history. Badoiu et. al. 1 gave an algorithm that computes the minimum 
enclosing ball approximately, with radius at most 1 + e times the optimal radius in time 0(nd/e 2 + 
(1/e) 10 ), independent of the number of dimensions, using convex programming. Their algorithm 
was based on the idea of Core- Sets, a small set of poly '(1/e) points whose minimum enclosing ball is 
almost as large as that of all the n points. This was improved to 0(nd/e + (1/e) 5 ) in [3| by finding 
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smaller core-sets of size |~l/e~|. They also provide a simple 0(nd/e 2 ) time algorithm for finding the 
minimum enclosing ball that does not require convex programming. Combining the two results 
gives a 0(nd/e + (1/e) 5 ) time algorithm for finding the minimum enclosing ball while eliminating 
the use of convex programming. 

For the minimum enclosing cylinder problem, Har-Peled and Varadarajan [0] gave an algorithm 

Of 1 log 2 - ) 

with running time of 2 v ^ c 'nd that finds a cylinder with radius at most 1 + e times the optimal 
radius. They also generalized their algorithm to computing the minimum radius /c-dimensional flat 

0(fc 2 ) 

in time exp( e ik+3 )nd, where the radius of a /c-flat is the maximum distance of the given set of 
points from this /c-flat. 

In this paper, we present algorithms for computing the minimum enclosing polytope for a large 
family of shapes, provided either only translations and no rotations are allowed, or only rotation 
about a fixed point is allowed; that is, one is allowed to only scale and translate a given shape, 
or scale and rotate the shape around a fixed point. We hope that it may be possible to combine 
the techniques for translation and rotation to solve the problem without these restrictions. Our 
algorithms are based on a simple greedy principle applied iteratively: simply move the current 
polytope directly towards any outside point till it touches the surface. 

For computing the minimum enclosing ball, this gives a simple greedy algorithm that repeatedly 
moves a ball directly towards the farthest uncovered point till it touches the surface. If we start with 
a ball of the optimal radius, we show that running 0(1/ e) such steps gives the optimal position of the 
ball approximately, within a running time of 0(nd/e) (section |2j) . This simple principle generalizes 
to arbitrary convex shape when only translations are allowed, requiring at most 0(1/ e 2 ) iterations 
(section EJ. It also works if the shape can be expressed as a union of a few convex shapes - however, 
requiring a running time exponential in 1/e 2 . Our algorithm implies that core sets of size 0(l/e 2 ) 
exist not only for the minimum enclosing ball but also for any convex shape with a fixed orientation. 

Next we look at covering a set of points by a convex body while allowing only rotation about 
a fixed point (section 0J). For a certain class of convex bodies with an axis of symmetry that 
includes cylinders, cones and ellipsoids, we prove that our techniques work provided the problem is 
confined to a half space bordering at the point of rotation. Without this restriction, we obtain an 
algorithm whose running time is exponential in 1/e 2 . This gives us an 2°( 1//e ^nd time algorithm 
for the min-cylinder problem provided we are given a fixed point on the axis. This also implies that 
core-sets whose size depend only on e exist for rotational problems as well. Although we are unable 
to combine our techniques for translations and rotations for general shapes, for the min-cylinder 
problem, we give an algorithm almost identical to the one in but with an improved running 
time of 2°^^° s ~^nd (section |SJ) . This generalizes to computing the minimum radius /c-dimensional 

n / p°( fc2 ) 1 \ 

fiat m time exp( — ^ — log - )nd. 

2 Minimum Enclosing Ball 

Given a set S of n points in d dimensions, we provide an algorithm to compute the minimum 
enclosing ball with radius at most 1 + e times the optimal in time 0(nd/e). 

We start with a simple algorithm MEB (figure ^) that works as follows: the algorithm starts 
with an arbitrary ball of the optimal radius, and for any point at least e outside this ball, moves 
the ball till the surface touches the outside point. This involves guessing the optimal radius of the 
minimum enclosing ball. Assume without loss of generality that the optimal ball is of unit radius. 

Let d(P,Q) denote the distance between two points P and Q. Let B(C,r) denote the ball of 
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Algorithm MEB 

1 Start with a ball of optimal radius. 

2 Repeat until every point is within 1 + e of the current center C. 

3 Find the farthest point P from C. 

4 Move C towards the point P till P touches 
the unit sphere centered at C. 



Figure 1: A simple algorithm for finding the minimum enclosing ball 




Figure 2: The center of the ball gets closer to the optimal position in each iteration 



radius r centered at point C. 

Theorem 1 Algorithm MEB terminates in 0(1/ e) iterations. 

The basic idea behind this theorem is that in each iteration the center C in the algorithm moves 
closer to the optimal center Copt- If <k is the distance of C from Copt in the i th iteration, then 
di decreases as follows. 

4 + x <d 2 -e 2 

This is because LCi,Ci + \,CoPT is obtuse (figure HJ), as PCi + \ is not shorter than PCopt 
implying lPCi + \CoPT is acute. So, df +1 < df — d(Ci, Cj+i) 2 <d\ — e 2 . Since the initial value of df 
is at most a constant, and since it decreases by e 2 in each iteration, the algorithm must terminate 
in 0(l/e 2 ) iterations. A tighter analysis based on the following lemma will show that it actually 
terminates in 0(1/ e) iterations. 

Lemma 2.1 If C and Copt o/re distance d apart, there must be a point that is at least d 2 /A from 
the surface of the unit ball centered at C. 
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Proof Look at the hemisphere in the optimal ball that is directly facing away from C. It is well 
known (for reference see lemma 5.2 in 4 ) that one of the n points, P, must lie on this hemisphere. 
Since LCCoptP is obtuse, 

d(C, P) > \/l + d 2 > 1 + d 2 /A ( since d < 2) 

This means that the point P is at least d 2 /4 away from the surface of the unit ball centered at 
C. m 

Now that we have a lower bound of d 2 /4 on d(Ci, Cj+i), we are ready to prove that algorithm 
MEB terminates in 0(l/e) iterations. 

Proof of theorem [Q Note that d 2 i+l < d 2 - d(d, C l+1 ) 2 < d 2 - df/W 

Let $j = df, and we get the recurrence relation <3?i+i = — <3? 2 /16. It is easy to check that 
if we start with $o = 1, then after 0{2 l ) iterations, <3?j decreases to < 1/2* (if <3? is 2~*, in 0(2*) 
iterations it will become less than 2 - (* +1 ) ). So after 0(1/ e) iteration, <3?j becomes < e. After that, 
since it decreases by at least e 2 in each iteration, there can be at most e/e 2 = 1/e further iterations. 
■ 

Algorithm MEB requires guessing the optimal radius. The distance of the farthest point from 
any given point in the set is within factor 2 of the optimal radius. A binary search with at most 
0(log(l/e)) tries can be used to ascertain the correct radius approximately within a factor of 1 + e. 
If a guess is too small the algorithm will not terminate in 0(1/ e) iterations. For a certain guess, if 
the algorithm terminates successfully, this means the guess is greater than a 1 + e approximation of 
the optimal radius and so the guess may be decreased. So by running MEB at most 0(log(l/e)) 
times the minimum enclosing ball can be computed in time 0(^p log |). 

2.1 Eliminating the Binary Search 

To eliminate the binary search, in algorithm MEBOPT (figure we start with a ball of radius 
less than optimal and increase it in certain iterations. The radius of the ball is always less than the 
optimal radius but the gap decreases as the iterations proceed. 

The algorithm maintains a lower bound r, and an error 5, such that r < roPT < r + 5. Again as 
before we assume that the optimal radius is 1. In each iteration we reduce the error 5 by a factor 
of 3/4 based on the following lemma. 

Lemma 2.2 After 0(1/ 5) iterations of moving the ball to farthest outside point, every outside 
point must be within 35 from the surface. 

Proof Let hi denote the distance of farthest point from surface of B(Ci,ropr), after i iterations. 
Then the farthest point is at a distance hi — 5 from B(Ci,roPT)- Just as in proof of theorem^ it 
is easy to check that (only difference is that d(Ci, Cj+i) = hi — 5 ) 

d 2 +i<d 2 -(^-<5) 2 

Now as long as hi > 25, we have df +1 < df — /i 2 /4. Again lemma l2~T1 savs hi > d 2 /A, implying 
that as long as hi > 25, we have d 2 +1 <d 2 — c^/64 So as in proof of theorem^ in 0(1/5) iterations, 
either di < 5 or hi < 25. In either case, distance of the farthest point from the center is at most 
r opt + 25 < r + 35. 
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Algorithm MEBOPT 


l 


Initialize C = any arbitrary point of the n points 


2 


r = 1/2 of distance of farthest point from any one point 


Q 
O 


5 = 1/2 of distance of farthest point from any one point 


A 

4; 


rvepeai until o ^ e 





For 0(1/5) iterations 


6 


find farthest point P from C 


7 


move B(C,r) till its surface touches P. 


8 


Let s = distance of the farthest point from the surface of current ball B 


9 


If s < 35/4 


10 


5 = 35/4 


11 


else r = r + 5/4 


12 


5 = 35/4 



Figure 3: An algorithm that does not require binary search 



Using the above lemma, we can test if the current estimate r in fact has an error of at most 5/4 
in 0(4/ 5) iterations. If after so many iterations, the farthest point is more than 35/4 away from 
the surface, then we can conclude that the error was more than 5/4 and so r can be increased to 
r + 5/4. Otherwise, since every point is within distance 35/4 outside the surface, we can conclude 
that the error in r is at most 35/4. Since each iteration in step 4 runs in time 0(nd/5), and 5 
decreases geometrically to e, we have proved the following theorem. 

Theorem 2 Algorithm MEBOPT finds an approximate minimum enclosing ball in time 0(nd/e) 

3 Generalizing to Convex Polytopes 

The simple algorithm of moving towards the outside point works not only for finding the minimum 
enclosing ball but also for minimum enclosing polytope of any given convex shape with a fixed 
orientation. That is, one is only allowed to translate and scale the given shape but is not allowed 
to rotate it. Again, for ease of exposition, we assume that the maximum inter-point distance is at 
most 1. We present algorithm MINCON (figure similar to MEB, that finds an approximate 
optimal solution in 1/e 2 iterations. Again as in algorithm MEB, we guess the optimal size of the 
given shape but do not know its position to begin with. We repeatedly find an outside point at 
least e away from the surface and move the current polytope by the shortest distance till the point 
touches the surface. 

Theorem 3 Algorithm MINCON terminates in 1/e 2 iterations. That is, after so many iterations 
no point will be more than e outside the surface. 

Proof The proof is very similar to that of theorem 1. Let Q be the point on the current polytope 
closest to P (figure |SJ). Also, we know that P is in the optimal polytope. Let P' denote the 
corresponding point in the current polytope. That is, the vector P'P is the displacement of the 
optimal polytope from the current polytope. We will prove that IP'QP is obtuse. If not, there 
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Algorithm MINCON 

1 Start with a polytope guessed to be of optimal size positioned anywhere. 

2 Repeat until done 

3 Find any point P that is at least e away from the 
surface of the current polytope 

4 Find the point Q on the polytope that is closest to P. 

5 Move the polytope so that Q coincides with P. 

Figure 4: Algorithm for finding Minimum Enclosing Convex Polytope 




Figure 5: The polytope keeps moving closer to the optimal position in each iteration 
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is a point on the segment P'Q that is closer to P than Q. And since both P' and Q are in 
the current polytope, that point must also be within the current polytope. So Q cannot be the 
closest point to P, which is a contradiction. If d, denotes the distance of the current polytope of 
the optimal one in this ith iteration, then di = d(P',P). After the displacement by the vector 
QP, the polytope will be off from the optimal position by the vector P'P — QP = P'Q. So 
d 2 +1 = d(P',Q) 2 < d(P',P) 2 - d(Q,P) 2 < d 2 - e 2 This proves that the algorithm terminates in 
1/e 2 iterations. ■ 

This technique also extends to shapes that are a union of a small number of convex-shapes 

Theorem 4 Given a shape and orientation that can be expressed as a union of c convex-shapes, the 
smallest enclosing polytope with that shape and orientation can be computed within e approximation 
in time (P^ e ^nd 

Proof The algorithm is identical to MINCON , except that at each step we guess one of the 
convex bodies that contains the outside point and move the polytope till that convex body touches 
the point. ■ 

This in fact proves that Core-Sets exist not only for minimum enclosing ball but also for any 
convex shape with a given orientation. 

Definition 1 Given a set of points, S, and convex shape and orientation, we say that a subset T 
of S forms a Core-Set if the minimum enclosing polytope of T has every point of S within distance 
at most e outside its surface. 

Theorem 5 For a set of points with maximum inter-point distance 1, and for a given convex shape 
and orientation, there is a Core-Set of size 0(l/e 2 ). 

Proof Instead of starting with a polytope of the optimal size, we start with one just small enough 
so that it can never be positioned to have every outside point within distance e from its surface. 
Now we know that if we start with this size, then algorithm MINCON would never terminate in 
the 0(l/e 2 ) iterations it otherwise would have. We let the algorithm run for one more than 0(l/e 2 ) 
iterations and look at the 1 + 0(l/e 2 ) points that are visited. We will prove that these points form 
the required Core-Set. The minimum enclosing polytope of these l + 0(l/e 2 ) points must be larger 
than the one we started with as otherwise, by theorem |21 algorithm MINCON would not require 
more than 1 + 0(l/e 2 ) iterations on these points. Since the initial polytope can be chosen so that 
every outside point is within distance arbitrarily close to e from the surface of this initial polytope, 
we have proved the theorem. ■ 



4 Allowing Rotations 

So far we did not allow the convex polytope to be rotated and only allowed translations. In this 
section we prove that our techniques work if only rotation about a fixed point and no translations 
are allowed, provided certain conditions are met. 

Given a polytope that has an axis of symmetry (that is, every cross section along the axis is 
hyper-sphere of dimension d — 1) with the axis passing through the origin, and a set of points S, 
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Algorithm MINROT 

1 Start with the axis as any ray in the given half-space shooting from the origin. 

2 Iteratively find any outside point and rotate the axis by the smallest 
angle till the surface of the polytope touches the outside point. 

We assume that the distance between the outside point and the point 
on the surface it touches is at least e, as otherwise we are done. 

Figure 6: Algorithm for rotational problem in a half-space 

our goal is to rotate the axis till the points in S, are covered by the polytope. We will also assume 
that the following conditions are satisfied 

• All these points and the optimal polytope lie in a half-space with the bounding hyper-plane 
passing through the origin. 

• Any d-dimensional hyper-sphere centered at the origin intersects the polytope in a single 
hyper-sphere of dimension d — 1. This d — 1-dimensional hyper-sphere divides the original 
d-dimensional hyper-sphere into two disjoint regions. We also assume that the interior of the 
polytope intersects the <i-dimensional hyper-sphere in the smaller of these two regions, (this 
is similar to the convexity requirement in algorithm MINCON . In three dimensions this 
would mean that every sphere passing through the origin cuts the polytope in at most one 
circle. Also the interior of the polytope intersects the sphere in the smaller of the two regions 
on the sphere formed by the circle). 

Examples that satisfy these conditions are cylinders, half-cones, ellipsoids lying in a half-space 
with axis passing through the origin. Again, for ease of exposition, we assume that all points are 
at most at unit distance from the origin. Our algorithm MINROT (figure EJ) repeatedly rotates 
the polytope by the smallest angle so as to touch an uncovered point at least e outside the surface. 

Theorem 6 Algorithm MINROT terminates in 1/e 2 iterations. That is, after so many iterations 
no point will be more than e outside the surface. 

Proof The proof is very similar to that of theorem ^ Let 6 be the angle between the current axis 
and the optimal one. We will argue that the distance between corresponding points on the two axes 
on the units sphere centered at the origin decreases in each iteration. This distance d = 2sin(#/2). 

Let P be the outside point chosen in a certain iteration and Q be the closest point on the 
current polytope in terms of rotation required to move Q to P. Look at the sphere centered at the 
origin passing through P (figure [7J). This sphere intersects the current polytope in a hyper-circle, 
C, passing through Q. The point P lies within the optimal polytope. Let P' be the corresponding 
point in the current polytope. P' must be on the sphere inside the hyper-circular region C. Look at 
the great circle on the sphere passing through P and Q. Project all points onto the two dimensional 
space containing this great circle. Under this projection the hyper-circle C will become a segment 
QQ'. Since the points P, Q, Q' lie on a half circle IPQQ' is obtuse. Since, under the projection, 
P' lies in the minor segment formed by QQ f , LPQP' is also obtuse. Scale the distances so that the 
great circle is of unit radius. So df +1 = d(P', Q) 2 < d(P' , P) 2 - d{Q, P) 2 < dj - e 2 ■ 

Algorithm MINROT assumes that the polytope lies in a given half-space. We now provide an 
alternate algorithm for polytopes with an axis of symmetry passing through the origin without 
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Figure 7: The axis keeps rotating to the optimal position in each iteration 



Algorithm FULLROT 

1 Start with the axis as any ray from the origin. 

2 For each outside point P look at the sphere centered at the origin. 

3 The sphere intersects the polytope in at most two circles, 
C\ and C2 on different sides of the origin. 

4 Now the axis could be rotated to either touch C\ or C2 to P. 
Guess one of them and rotate the axis by the smallest 
angle till P touches the chosen circle. 

Figure 8: Algorithm without the half-space assumption 



the half-space assumptions in MINROT. We will assume that the polytope is symmetric around 
the origin and its intersection with any hyper-sphere is at most two equal sized disjoint d — 1- 
dimensional hyper-spheres on opposite sides of the origin. Again the interior of the polytope 
intersects the original hyper-sphere in the smaller of the two regions formed by each of the two 
d — 1-dimensional intersection hyper-spheres. Examples are cylinders, cones, ellipsoids centered at 
the origin. However this algorithm runs in time 2°( 1 / e2 W. For ease of exposition, this algorithm, 
FULLROT (figure is described for three dimensions. 

Theorem 7 Algorithm FULLROT terminates inl/e 2 iterations. That is, after so many iterations 
no point will be more than e outside the surface. The deterministic version of this algorithm runs 
in time 2 0<yl / € ^nd by trying all possible guesses. 

Proof As in the proof of theorem 4 in any iteration there would be two points Q\ and Q2 on 
the circles C\ and C2 closest to the outside point P. As before P' would lie in one of the minor 
segment formed by one of C\ and C^- We guess the correct one, say C\. Again project all points 
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to the plane containing the great circle passing through P and Q\. Since the angle Q2Q1Q1 is 90 
deg, the angle PQiQ[ is obtuse. The rest of the proof is same as that of theorem El ■ 

Again, as before, we can extend our techniques to shapes that are a union of a small number 
of bodies that satisfy the conditions required by algorithm FULLROT. 

Theorem 8 Given a shape and orientation that can be expressed as a union of c shapes, each 
satisfying the conditions required by algorithm FULLROT, we can find the smallest enclosing 
polytope with that shape and orientation in time (2c)°^ 1 ' e ^nd 

Just as in section |21 we can derive core-sets for rotational problems. 

Theorem 9 For rotational problems with shapes satisfying conditions for algorithms MINROT 
and FULLROT, core-sets of sizes 0(1/ e 2 ) and 2°( 1 / e ' exist, respectively. 

Note that an infinite cylinder with its axis passing through the origin satisfies the assumptions 
of algorithm FULLROT. So we have: 

Corollary 1 For the minimum radius cylinder problem if we are given a point on the axis of the 
optimal cylinder, algorithm FULLROT runs in time 

Note that in the min-cylinder problem, the maximum distance between all points may not be 1 
as assumed. This can be easily overcome by setting the initial position of the axis to pass through 
the farthest point from the origin - we omit the details here. 

5 Minimum Radius Cylinder 

Although we do not have any general results for a combination of rotation and translation for 
different shapes, we provide an algorithm for the min-cylinder problem without restrictions that 
runs in time 2°( 1//e ^nd. The algorithm is similar to the one mentioned in [S] with a running 

Of 1 log 2 -) 

time of 2 ^ e 'nd. Our algorithm can be viewed as following the greedy principle underlying 
the other algorithms of this paper: In each iteration it moves the axis of the cylinder along the 
plane containing the axis and an outside point by guessing its optimal position in that plane 
approximately. 

Without loss of generality assume that the optimal radius is 1. We will show later how this 
optimal radius can be computed approximately using a binary search. We start with a certain 
initial position of the axis that will be specified latter. Let Iopt be the axis of the optimal cylin- 
der. Let U and V be the farthest two points on Iopt that are projections of points in set S on Iopt- 

Algorithm MINCYN - 

1. We iteratively compute an estimate I4 of Iopt and points f7, and Vi on l{ that are close to 
projections of U and F on ij. 

2. We maintain the following invariant: d(U{, U) < 5 and d(Vi, V) < 5 

3. Zj+i , J7i+i and Vi+i are computed from l%,Ui and Vi as follows: Find any point P that is at 
distance more than 1 + e from Zj. Look at the plane h containing li and P. We will try to 
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set Ui+i and V, + ± close to Uh = proj(U, h) and Vh = proj(V, h) respectively, by the following 
process. 

From the invariant, we have d(Ui, Uh) < 5. So Uh lies in a circle in h of radius 5 centered at 
Create a mesh, where each element has side e/8, so that U% itself is a mesh point, and 
guess the mesh point closest to Uh and at a distance at most 5 from U\. We need to guess 
one out of j^rrbi points and set this point to Ui + \. Clearly this point is at most e/4 from Uh- 
Similarly we guess V,+i out of at most 0(1/ e 2 ) points. 

We will prove convergence by arguing that the potential function, $ = d(U,Ui) 2 + d(V,Vi) 2 , 
decreases significantly in each iteration. 

Lemma 5.1 We maintain the invariant, d(Ui,U) < 5 and d(Vi, V) < 5, during each iteration of 
algorithm MINCYN. 

Proof We will show that d(Ui, U) only keeps decreasing and if the invariant is true to start with, it 
always remains true. Now for any point X on the plane h, d(X, U) 2 = d(X, Uh) 2 + d(Uh, U) 2 . Since 
we choose Ui+i to be the mesh point closest to Uh, among mesh points including f/j, d(Ui+\, Uh) < 
d(Ui,Uh)- So the invariant 2 follows. ■ 



Lemma 5.2 In each iteration of algorithm MINCYN , the potential function $ decreases by at 
least e 2 /2 

Proof Note that UUh is perpendicular to the plane containing Uh,Ui and JJj+i. So, 

d(U,U i+1 ) 2 = d(U,U h ) 2 + d(U h ,U l+1 ) 2 

= d(U,U t ) 2 -d(Ui,U h ) 2 + d(U h ,U i+1 ) 2 
< d(U,Ui) 2 -d(Ui,Uh) 2 + e 2 /A 

Similar inequality holds for d(V, Vi + i) 2 . Adding the two we get, < 3>i — d(Ui,Uh) 2 — 

d(V h V h ) 2 + e 2 /2. 

We will prove that at least one of d(Ui, Uh) and d(Vi, Vh) is more than e. For if not then we will 
show that P cannot be within distance 1 of any point in the segment UhVh, which is a contradiction 
because UhVh is the projection of UV. 

Let l p be the line in plane h passing through P and perpendicular to Zj, meeting Zj at P' . Then, 
d(P, P') > 1 + e. Project all points to l p . The segment UhVh projects down to a segment U p V p . 
Since there is a point on UhVh that is at most distance 1 from P there must also be such a point, 
Q, on UpVp. Since d(P', Q) > e, at least one of U p and V p must be at least e away from P' . Since 
distances only decrease under projections, at least one of d(Ui,Uh) and d(Vi,Vh) must be > e. 

So, we get < $j - e 2 + e 2 /4 < $ 4 - e 2 /2 m 

Finally we need to prove that we can choose an initial line Iq and points on it Uq and Vb that 
satisfy the invariant. Look at any point X in the set S. Let Y be the farthest point from X in S. 
Set Iq to the line passing through X and Y. It is easy to verify that every point in S must be within 
distance 4 from Iq (See lemma 5.2 in jU] for a more general statement). Look at the projections of 
points in S on Iq and I. For any point Z in S, let Zq denote its projection on Iq and Zi denote its 
projection on I. Then d(Z, Zq) < 4 and d(Z, Z{) < 1. So d(Zo, Z{) < 5. 
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Set Uq, and Vo to the farthest two points among projections of points of S on Iq. Clearly these 
points are at most at distance 5 from U and V respectively. 
So we have proved the following theorem 

Theorem 10 Algorithm MINCYN terminates in 0(1/ e 2 ) iterations. That is, after so many 
iterations no point will be more than e outside the surface of the cylinder. A deterministic version 
of this algorithm runs in time 2°^ log ~^nd. 

The deterministic version follows by simply eliminating the guess. Each guess requires guessing 
twice from 0(1/ e 2 ) choices, this guessing happens at most 0(1/ e 2 ) times. We also need to clarify 
how the optimal radius required by the algorithm can be determined. The distance of the farthest 
point from the initial position of the axis Iq is within a constant factor of the optimal radius. The 
algorithm MINCYN terminates only if the radius used is larger than the optimal radius and does 
not terminate if the radius used is too small. A binary search involving log(l/e) trials will result 
in a value that is within 1 + e of the optimal radius. 

Using the techniques in [§] this algorithm generalizes to computing the min-radius /c-dimensional 
flat - we omit the details here. 

Theorem 11 The minimum radius k- dimensional flat can be computed in time 2 ^ e nd. 

Acknowledgments 

I would like to thank Mihai Badoiu, Piotr Indyk, Sariel Har-Peled and Kasturi Varadharajan for 
useful discussions. 

References 

[1] Mihai Badoiu, Sariel Har-Peled and Piotr Indyk. Approximate clustering via core- sets. Proceed- 
ings of the 34th Symposium on Theory of Computing, 2002: 250-257. 

[2] M. Badoiu and K. L. Clarkson. Optimal core-sets for balls. Manuscript. 

[3] Mihai Badoiu and Ken Clarkson. Smaller Core-Sets for Balls. Proceedings of the 14th Annual 
ACM-SIAM Symposium on Discrete Algorithms (SODA'03): 801-802. 

[4] Ashish Goel, Piotr Indyk, and Kasturi R. Varadarajan. Reductions among high dimensional 
proximity problems. Proceedings of the 12th ACM-SIAM Symposium in Discrete Algorithm, 
2001: 769-778. 

[5] P. Gritzmann and V. Klee. Computational complexity of inner and outer j-radii of polytopes in 
finite- dimensional normed spaces. Math. Program., 59:163213, 1993. 

[6] P. Gritzmann and V. Klee. On the complexity of some basic problems in com-putational con- 
vexity: I. containment problems. Discrete Math., 136:129174, 1994. 

[7] M. Grotschel, L. Lovasz, and A. Schrijver. Geometric Algorithms and Combi-natorial Opti- 
mization, volume 2 of Algorithms and Combinatorics. Springer- Ver lag, Berlin Heidelberg, 2nd 
edition, 1988. 2nd edition 1994. 



12 



[8] S. Har-Peled and K. R. Varadarajan. Projective clustering in high dimensions using core-sets. 
In Proc. 18th Annu. ACM Sympos. Comput. Geom., pages 312318, 2002. 

[9] Sariel Har-Peled, and Kasturi R. Varadarajan. High-Dimensional Shape Fitting in Linear Time. 
Proceedings of the Symposium on Computational Geometry, 2003: 39-47. 



13 



